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We use recent lattice data on the heavy quark potential in order to determine the dissociation 
temperatures of different quarkonium states in hot strongly interacting matter. Our analysis shows 
in particular that certain quarkonium states dissociate below the deconfinement point. 



I. INTRODUCTION 



The behavior of the heavy quarkonium states in hot strongly interacting matter was proposed as test of its confine- 
ment status, since a sufficiently hot deconfined medium will dissolve any binding between the quark- antiquark pair 
0. Another possibility of dissociation of certain quarkonium states (subthreshold states at T = 0) is the decay into 
open charm (beauty) mesons due to in- medium modification of quarkonia and heavy- light meson masses Q. 

The production of J and T mesons in hadronic reactions occurs in part through production of higher excited cc 
(or hh) states and their decay into quarkonia ground state. Since the lifetime of different subthreshold quarkonium 
states is much larger than the typical life-time of the medium which may be produced in nucleus-nucleus collisions 
their decay occurs almost completely outside the produced medium. This means that the produced medium can be 
probed not only by the ground state quarkonium but also by different excited quarkonium states. Since different 
quarkonium states have different sizes ( binding energies ), one expects that higher excited states will dissolve at 
smaller temperature than the smaller and more tightly bound ground states. These facts may lead to a sequential 
suppression pattern in J /ip and T yield in nucleus-nucleus collision as the function of the energy density. 

Here we will discuss heavy quarkonium dissociation below the deconfinement point where it is due to in-medium 
modification of the open charm (beauty) threshold 1^ as well as above the deconfinement point due to the well known 
screening phenomenon |^. 



II. QUARKONIUM PRODUCTION AND FEED-DOWN 

It is well known that J production in hadron-hadron collision is to a considerable extent due to the production 
and subsequent decay of higher excited cc states • The feed-down from higher excited states was systematically 
studied in proton- nucleon and pion-nuclcon interactions with iOOGeV incident proton (pion) beams In these 
studies the cross sections for direct production of different charmonium states (excluding feed-down) were measured. 
Then making use of the known branching ratios B[xi{lP) ip{lS)] = 0.27 ± 0.02, B[x2{lP) ^ ip{lS)] = 0.14 ± .01, 
and B[ip{2S) ip{lS)] — 0.55±0.05, one obtains the fractional feed-down contributions fi of the different charmonium 
states to the observed J/ip production; these are shown in the second and third columns of Tab. 1 

In the case of bottomonium the experiment provides only the inclusive (i.e. including also the feed-down from 
higher states) cross section for different (nS) states 1^ . The feed-down from (nP) states is known only for transverse 
momenta pr > 8GeV/c To analyze the complete feed-down pattern, we thus have to find a way to extrapolate 
these data to pr — as well as to determine the direct cross section for different {nS) states. This can be done using 
the most simple and general model for quarkonium production, the color evaporation model In particular this 
model predicts that the ratios of cross sections for production of different quarkonium states are energy independent. 
This prediction was found to be true for a considerable range of energies ]Tc[ |. The ratios between the different xii^P) 
states in this model are predicted to be governed essentially by the orbital angular momentum degeneracy Q; we 
thus expect for the corresponding cross-sections 

Xo(lP):Xi(lP):X2(lP) = l:3:5. (1) 

From Table 1 we have for tt^A^ collisions X2(l-P)/xi(l^) — 1-44 ± 0.38 and thus reasonable agreement with the 
predicted ratio 1.67. Actually, for pN interactions, the experiment measures only the combined effect of xi and X2 
decay (30 % of the overall J/ip production); the listed values in Tab. 1 are obtained by distributing this in the ratio 
3:5. 

Using considerations based on color evaporation model, in particular Eq. (1), the feed-down from higher excited 
bb states to T production can be predicted Q; the feed-down fraction are summarized in Tab.l. Alternatively the 
feed-down fraction from higher excited bb states can be predicted using NRQCD factorization formula Q. The results 
of this analysis are summarized in the last column of Tab. 1 
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TAB 1: Feed-down fractions for from higher excited states to the J jil) and T states. 



III. QUARKONIUM DISSOCIATION BELOW DECONFINEMENT 

Recent lattice calculations of the heavy quark potential show evidence for the string breaking at finite temperature 
1^. On the lattice the potential is calculated from the Polyakov loop correlator, to which it is related by 

y (r, T) = - In < L{r)L^ (0) > +C, (2) 

where L(r) is the Polyakov loop (see e.g. Ref. |^ for definition). The normalization constant C contains both the 
cut-off dependent self-energy and the entropy contributions ~TS (for T ^ — In < L(r)L'l'(0) > is actually the free 
energy of the static QQ pair). For a properly chosen normalization constant C, V(t,T) is the ground state energy 
of an infinitely heavy QQ pair. In absence of dynamical quarks (quenched QCD) V(r^T) is linearly rising with r for 
large separations indicating the existence of a flux tube (string). If dynamical quarks are present the flux tube can 
decay (the string can break) by creating a pair of light quarks qq from the vacuum once V(r^ T) is larger then twice 
the binding energy of a heavy-light Qq (qQ) meson Q Thus the potential at very large distances is constant VooiT) 
and is equal to twice the binding energy of a heavy-light meson. 

It is expected that medium effects are not important at very short distances. Therefore at very short distances the 
potential V{r, T) should be given by the Cornell potential ||l^ 

V(r) = --+ar (3) 
r 

We use this fact to determine the normalization constant C and set the potential to be of the Cornell form at the 
the smallest distance rT = 0.25 available in lattice studies of with e = 0.4 as expected for (2-|-l)-flavor QCD p2[ . 
The resulting potential and Voc{T) are shown in Fig. |l|. Note the strong temperature dependence of Voc{T). Since 
for sufficiently heavy quarks {rriQ 3> ^qcd) it does not matter whether the quark is inflnitely heavy or just merely 
heavy, the open charm (beauty) meson masses are approximately given by 2A/£).s(r) = 2mc,b + Voo{T). 

Now the temperature dependence of the different quarkonium states should be addressed. At zero temperature the 
heavy quark masses permit an application of potential theory for description quarkonium spectroscopy (see e.g. |p4|). 
Furthermore it turns out that the time scale of gluodynamics relevant for quarkonium spectroscopy is smaller than 
(mgw) ^ (v being the heavy quark velocity) ||Tl| ]. For sufficiently heavy quarks this time scale is much larger than the 
typical hadronic time scale A^^^ ^ l/"^. The decay of the flux-tube like all other hadronic decays has time scale of 



^Similar phenomenon occurs of course at T = 0. However it is much more difficuh to observe it on lattice (see e.g. [[Ll[). 
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order 1/m. Therefore in the potential theory the potential must always be of Cornell form (i.e. linearly rising at large 
distances). These considerations have direct phenomenological support. Namely, simple potential models with linearly 
rising potential can describe reasonably well also the quarkonium states above the open charm (beauty) threshold. 
Many of this higher excited states have effective radius of order or even larger than Ifm [[l2yi3[ . Contrary to this 
situation in the case of the potential becoming flat around 1/m (the expected radius of string breaking at T = 0) the 
higher excited states above the open charm (beauty) threshold simply do not exist. Therefore we have determined the 
temperature dependent heavy quarkonium masses from the Schrodingcr equation with the temperature dependent 
string potential 
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FIG. 1. The heavy quark potential and its asymptotic value below deconfinement at different temperatures. The line on the 
left figure is a fit to the data points. 



Vstring{r,T) = -(e - iarctan(2rr)) ^ + 

(a(r) - ^ - ^arctan-^) r + ^ ln(l + Ar^T^) 



(4) 



This form of the potential describes quite well the temperature dependence of the heavy quark potential in quenched 
QCD for appropriately chosen (t{T) In order to make contact to real QCD we set as = 0.4 and use Tcl\fo = 

0.425 from H for the deconfinement temperature ( by a we always denote the string tension at zero temperature). 
Furthermore we use the following values of the heavy quark masses, mc = l.SGeV and mi, — A.72GeV as well as 
•y/ff = OAAGeV for the zero temperature string tension. This values of parameters give a fairly good description of 
the observed quarkonium spectrum at zero temperature. The temperature dependence of the string tension was taken 
from [Tsl ] . The resulting quarkonia masses are shown in Fig ^j. 

Since the smallest distance available on lattice is only 0.25T^^ one may worry about possible medium effects at 
this distance and their role in determination of 14o(T). Normalizing the Polyakov loop correlator (H) at r = 1/(4T) 
with Eq. (U) , we thus obtain what might be a more reliable estimate of the plateau Voo {T) than with the T = form 
(^. It turns out, however, that the two forms of short distance behavior resulting from the zero temperature Cornell 
potential (H) and (jj) are practically identical, so that the normalization is in fact not affected by the in-medium 
modifications at larger distances. To consider further possible uncertainties of the normalization procedure, we have 
also normalized the Polyakov loop correlator at the next smallest distance r — \/2/(4T). The resulting two forms 
of Voo{T) are shown in both Fig. |2| The difference between the two curves of Vcxi{T) provides an estimate of the 
normalization error. Except for the region very near T = Tc, the uncertainty is seen to be quite small. 

From Fig. || one can see that ip' and Xc states become an open charm states well below and can dissociate by 
decaying into DD. The situation is similar for T(3S') and Xbi'^P) states which can decay into BB below Tc- For J/tp, 
Xbi^P) and T(25') it is not possible to say whether they will dissociate above Tc or just below Tc- Finally, the T(IS') 
state will definitely dissociate above the deconfinement. 



IV. QUARKONIUM DISSOCIATION BY COLOR SCREENING 

In the deconfined phase it is customary to choose the constant C in (||) to be the value of the correlator at infinite 
separation C = Tin < L{r)L^{0) >= \ < L > \^. The resulting connected correlator defines the so-called color 
averaged potential 

n.,r)^-Tlnl^i!:l^ (5) 
I <T> I' 
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FIG. 2. The masses of different quarkonia states and the open charm (beauty) threshold as function of the temperature. 
Shown are the charmonia masses and open charm threshold (left) and bottomonia masses and open beauty threshold (right) 
as function of the temperature. The thick solid line is the open charm (beauty) threshold obtained from normalization at 
r = l/(4r). The thin solid line in the open charm (beauty) threshold obtained from normalization at r = \/2/{iT) (see text). 



The color average potential can be written as the thermal average of the potentials in color singlet Vi (T, r) and 
color octet Vs{T,r) states: 

exp(-T/(r,r)/T) = icxp(-^i(r,T)/T) + ^exp(-y8(r,T)/T) (6) 

In potential models it is assumed that quarkonium is dominantly a singlet QQ state. Furthermore the octet channel 
is repulsive (at least in perturbation theory) and therefore only a singlet QQ pair can be bound in the deconfined 
phase. Thus we need to know the singlet potential. The lattice data in the relevant case of 3 flavor QCD exist only 
for the averaged potential The averaged potential in 3 flavor QCD are shown in Fig. ^ for three representative 

temperatures. Note that within present accuracy of the lattice calculations the potential vanishes beyond some 
distance roiT) denoted by vertical arrows in Fig. ^ In perturbation theory, the leading terms for both are at high 
temperature and small r (r << T~^) of Coulombic form, 

V,{T,r) = --^, Vs{T,r)^+-^, (7) 

with a{T) for the temperature-dependent running coupling. In the region just above the deconflnement point T = Tc, 
there will certainly be significant non-perturbative effects of unknown form. We therefore first consider the high 
temperature regime, which we somewhat arbitrarily define as T > 1.45 Tc. In this region, we attempt to parameterize 
the existing non-perturbative effects through a conventional screening form, replacing Eq. (^ by 

- Vi(T,r) = 6 Vs{T,r) ^ ^ exp{- pi{T)r} , (8) 
4 r 



where fi{T) denotes the effective screening mass in the deconfined medium. We fit the lattice data by (g/|§) assuming 
a{T) and /i(T) to be unknown functions of T. Such a fit can describe the lattice data for T > 1.45Tc very well. 
Furthermore the temperature dependence of a{T) can be well described by 1-loop formula for the coupling in QCD 
with Aqcd = (0.34 ± 0.01)Tc and the screening mass /i(r) turns out to be constant in units of the temperature 
u(T) = (1.15 ± 0.02)T. A similar behavior of the screening mass was found in pure SU(2) and SU(3) gauge theory 
|l8| - . We shall now assume that the above form of the screening mass continues to remain valid as we lower the 
temperature to Tc- Such a constant screening mass down to Tc is again expected from studies of pure gauge theory 
[Q. On the other hand, quenched QCD (pure SU(3) gauge theory) studies indicate that when T is lowered to Tc, 
the perturbative ratio Vi/Vs = —8 will increase in favor of the singlet potential ||2l[]. We therefore try to describe the 
behavior just above Tc by a potential of the form (^), in which the color octet potential is given by 

Vs{T,r)^^^eM-f^r} (9) 
b r 

instead of Eq. (^; the factor c(T) < 1 accounts for the expected reduction of octet interactions as T ^ Tc- In the 
interval Tc < T < 1.45 Tc we thus fit the lattice results for V{T, r) in terms of the parameter c(T) with a{T) and /i(T) 
given above. In this way we get again quite a good fit of the lattice data Refs. P,0 temperatures T < 1.45Tc. 
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Now we are in a position to discuss quarkonium dissociation due to color screening. It is natural to assume that 
the heavy QQ pair cannot exist as a bound state if its effective binding radius ( the mean distance between Q and Q) 
is larger than the screening radius of the medium. The effective radii for different bound states are calculated from 
the Schrodinger equation with the singlet potential described above. 
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The screening radius of the medium can be identified with l/fi{T). However, the value of fi{T) strongly depends 
on assumption we have made to determine it. A more model independent and more conservative approach would be 
to identify the screening radius with ro(T) defined above. We use the latter approach. In Fig. ||we show the effective 
radius of J/ip and T states and the screening radius as function of the temperature. The intersection of these curves 
defines the dissociation temperature of J/ijj and T states. Similar analysis was done for excited states which may 
survive above Tc. 



V. SUMMARY AND CONCLUSIONS 



We have considered quarkonium dissociation in hot strongly interacting matter below as well as above the decon- 
finement. In a confined medium dissociation of certain quarkonium states occur due to in- medium modification of the 
open charm (beauty) threshold as well as the quarkonia masses. In the deconfined medium quarkonium dissociation 
is due to color screening. We summarize the dissociation temperature of different quarkonium states in Tab. 2. 
Combining these dissociation temperature with the feed-down fractions determined in section II we can predict the 
sequential suppression pattern of J/i/' and T states as function of the temperature. These are summarized in Fig. |^. 
For a more accurate determination of the quarkonium suppression patterns, it would be desirable to carry out direct 
lattice studies of the color singlet potential and of its quark mass dependence, which may become important near the 
critical temperature. Furthermore, to make contact with nuclear collision experiments, a more precise determination 
of the energy density via lattice simulations is clearly needed, as is a clarification of the role of a finite baryochemical 
potential. For the latter problem, lattice studies are so far very difficult; nevertheless, a recent new approach p^ ] 
could make such studies feasible. 

The methods used in the present study of in-medium modification of hadrons containing heavy quarks is of limited 
validity. A more accurate study of hadron properties in medium should be based on determination of hadron spectral 
function and their temperature dependence. Recent lattice studies indicate that a determination of hadron spectral 
function is possible at least within quenched approximation with the present days computer resources [ p3[ . Although at 
present the hadronic spectral function can be determined only in quenched approximation they may provide valuable 
test for the approach described here. In principle such studies will also include the dissociation of quarkonium states 
due to interaction with partonic constituents (thermal activation) in terms of width of the peaks in the calculated 
spectral function which is not accessible within present approach. 
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TAB 2: Dissociation temperatures of different quarkonium states. 
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FIG. 4. The suppression pattern of J/?/; (left) and T yield (right) as function of the temperature. 
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